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HEAT  TRANSFER  TO  LAMINAR  FLOW  IN  A  ROUND  TUBE  OR 
FLAT  CONDUIT  ---  THE  GRAETZ  PROBLEM  EXTENDED 


By.  tin  Sellars,  Myron  Tribus,  and  John  Klein 


NOIIENCLATURB 


Tli«  following  nooiane loture  ie  uaad  In  th«  papar: 

•  Coefficient  occurring  In  Evietlon  (10). 

-  Coefficient  occun-lng  In  Equation  (10). 

-  &lf  vldth  of  flat  duet,  ft. 

•  Coefficient  In  Equation  (3). 

a  Unit  heat  capacity  at  conetant  preaaura,  Btu/lb>*F. 

•  Coefficient  occurring  la  Equation  (IS). 

•  Coefficient  occurring  In  Equation  (18). 

•  F(a,  r*),  laplace  traaafora  of  Oraett  aolutlon. 

•  g(X*',  r*)  Integrating  kernel  for  baat-fluz  probleBi,  eee 
Equation  (Uk). 

•  0(a,  r'*')  Ijq>lace  tranefora  of  g. 

•  h(x^)  Integrating  kernel  for  heat-flux  prohlasB,  see 
Equation  (3U). 

•  H(s)  laplace  tranefom  of  h. 


Beseal  function  of  flrat  kind,  zero  order. 

Beaael  function  of  flret  kind,  l/3  order. 

Beaael  function  of  flrat  kind,  -l/j  order. 

Ibemal  conductivity  of  fluid,  Btu/aec-ft‘  (*r/ft). 
Coefficient  occurring  In  Equation  (A-3). 

Prandtl  aodulua,  dlaenalonleaa,  (uCp/k)(3fi00  g^). 

q(x)  heat  flux  per  unit  vail  area,  Btu/hr-ft*. 
laplace  tranefOm  of  (kq/rg). 


Radlua,  ft. 

Tube  radlua,  ft. 
(r/xo). 


SeyooUa  aodulua,  dlaenalonleaa,  (20^  rop/uge)  or  (W^bp/ugg). 

nranafora  variable. 

t(x'*’,  r*)  te^erature,  *T. 

T(b,  e*),  laplace  tranefora  of  t. 

Velocity  of  fluid,  ft/aec. 

Average  fluid  velocity  In  tube,  ft/aec. 

Ptataace  along  tube,  ft. 

(x/ro)(RePr)-i  or  (x/b)(HePr)-i. 


Dlatanee  froa  duct  vail,  ft. 


Dlatance  frca  tube  vail,  ft. 


Accesion  For 

NTIS  CRA&I 
DTIC  TAB 
Unannounced 
Justification 


ZeroofH(a).  . . . 

Elgenralue.  _ 

By _ 

viacosity  of  fluid,  ib-aec/fta.  Distribution  / 

Fluid  density,  lb/ft».  r-... 


Availability  Codes 


Tsrl«blea 
Sua^  TarlA^lea 


A-/ 


Avail  and/or 
Special 


IHTROIXJCTION 

The  problem  considered  here  is  posed  by  a  system  in  which  a  fluid  of 
constant  properties  flows  in  steady  laminar  motion  in  a  round  tube  or  flat  duct. 
The  velocity  profile  is  fully  established  and  peirabolic.  to  a  point  (x  «  0) 
the  fluid  is  isothermal.  After  this  point  a  prescribed  heat  flux  or  temperature 
is  given  at  the  wall  of  the  conduit  and  the  problem  is  to  find  the  temperature 
distribution,  as  well  as  the  connection  between  heat  flux  and  wall  tenperatvire . 
The  application  of  this  solution  to  practical  problems  of  heat  exchange  has  al¬ 
ready  been  so  well  established  that  further  comment  is  unnecessary. 


The  prohlem  has  been  considered  in  detail  by  a  nuniber  of  workers  and 
an  excellent  review  is  contained  in  the  book  "Heat  Treuisfer"  by  M.  Jakob.  ^  The 
problem  readily  reduces  to  the  finding  of  eigenvalues  €md  prior  to  this  paper 
only  the  first  three  eigenfunctions  and  the  first  four  eigenvalues  have  been 
known.  A  recent  papers  has  brought  out  the  lii?>ortance  of  obtaining  more  eigen¬ 
values,  euid  by  using  the  complete  set  of  eigenvalues  and  the  methods  of  refer¬ 
ence  2  the  classical  "Graetz  Problem"^  is  extended  to  more  cQ05)llcated  boundary 
conditions. 


The  problem  can  be  stated  in  mathematical  terms  as  follows: 


Given 


t  s  t 

'  ^  ^  A  >A.  ' 

M-  *  ^  [  I  ~  ] 


and  for 


with  either 


t-t. 


X  <  o 


or 


•t(x,  a.)  * 

^  A.)  s  ^(%) 


%  >o 


find  t(x,r)  and  the  relation  between  q(x)  and  t^(x). 

The  nondimens lonal  form  of  the  equations  is 

The  boundary  conditions  are 

t  •♦, 

Sind  either 

or 


X'^<  o 


■) 

V  %  1  J 


X^>o 


(1) 


(2) 


In  view  of  the  linearity  of  Equation  (2),  it  is  necessary  to  have 
the  fundamental  solution,  known  as  the  Graetz  solution,  to  construct  all  other 
needed  solutions.  Therefore,  the  Initlea  step  is  the  coopletion  of  the  Graetz 
solution. 

^  'Heat  Transfer,*  by  Max  Jakob,  John  Wiley  &  Sons,  Inc,,  New  York, 

N.  Y.,  vol.  1,  1949. 
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THE  CaiAETZ  SOnfTlON 


The  problem  considered  by  Graetz  and  most  other  vorkers  Is  Equation 
(2)  with  boundary  conditions. 

i  *1  x*<o 


%*  >o 


Led  0  be  a  solution  of  Equation  (2),  then 


e  =  Z 


.Kt* 


where  the  are  the  eigenvalues  required  to  make  the  solution  to  the  following 
differential  equation 


satisfying  the  boxindary  conditions  Rn(l)  •  0>  ^n(O)  *  coefficients 

are  determined  from  the  relation 


A*  (I- A*')  KU* 


-A 

/ 


^  S 


The  eigenfunctions  and  eigenvalues  have  been  given  only  for  n  =  1,  2, 
3.  The  higher  modes  of  Equation  (4)  are  very  difficult  to  calculate  for  large 
values  of  X  .  Therefore,  to  obtain  Xn  and  Cn  for  n  >  3,  a  solution  is  soiight 
which  will  be  valid  as  X^  ■*■  <».  It  will  be  found  that  the  resulting  formulae 
will  provide  good  answers  even  when  X^  is  small.  First,  look  for  a  solution  in 
the  form 

A  »  • 

and  find  that  gCr*^)  satisfies 

3"  ♦  j‘‘-t (6) 

Now  an  asynptotlc  solution  is  sought  in  the  form 

j  .  Xj.  .  j,  +>,-'3^....  (7) 

Substitution  in  Equation  (6)  and  equating  powers  of  X  gives 
3#  =  ±  t 


(8) 


4 


<J. 


(9) 


Since  \  is  large,  the  remaining  terms  in  Equation  (7)  are  neglected. 
Substitution  of  Equations  (8)  and  (9)  in  Equation  (7)  gives  for  R 


A  ®  ^36 

<7 


(10) 


Equation  (10)  is  the  so-called  WKB  approximation  auid  is  valid  for  o<r+  <  1  for 
sufficiently  large  Now  the  coefficients  A  and  B  mist  be  determined  so  that 

Equation  (10)  will  correspond  to  the  regular  solution  of  Equation  (4),  \diere  r+ 
is  small.  For  small  r"*"  Equation  (10)  is 

-iX»* 

»  .  A*  (U) 

(Hf  (I  -A-'*)*' 

Inspection  of  Equation  (4)  shows- that  when  r^  is  small  enough  so  that 
X.2  (1  _  r+®)  the  classical  solution  behaves  as  J©  (Ar+)#  since  Equation 

(4)  then  becomes  a  Bessel  equation.  For  large  Xr"*",  even  if  r"*"  is  small,  the  as- 
yn5)totic  expression  for  Jo(Xr'*‘)  is 

(12) 

and  thus,  it  is  seen  that  to  make  Equations  (11)  and  (12)  equal  for  r+  small,  it 
is  required  that 


8uad  for  0  <  r"*"  <  1 


Equation  (l4)  is  not  a  good  approximation  to  the  solution  as  r"*"  -►  1, 
since  it  has  a  singularity  there.  Because  a  boundary  condition  is  to  be  Imposed 
at  r+  =  1,  the  development  of  an  altenaate  solution,  valid  neeu:  r"*"  =  1,  is  con¬ 
sidered.  By  patching  it  on  to  Equation  (l4)  the  solution  over  the  range  0  A  r"*" 

^  1  is  obtained. 

The  following  change  of  variable  is  made 

I*  =  l-A* 

and  Equation  (4)  becomes 
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A  .  _I_  f  X  (I5; 

Tp  '-J' 

Now  consider  0  <  z"*"  «  1  and  define  a  new  variable 

H  -  X*’*  (16 

Substitution  of  Equation  (16)  into  Equation  (lU)  yields  for  large  X 

^  (17) 

which  has  the  solution 

A  •!>  ^ (l®> 

Hie  constants  D  and  E  are  to  be  so  chosen  that  for  small  s'**  Equations  (I8)  and 
(llv)  are  equivalent. 


Change  the  variable  from  r+  to  z+  in  Equation  (ill-)  and  perform  the  in¬ 


tegration 


( f— —  r '  ^ 

3.  M  ■  j.  /TT’ ♦  )■ 


For  small  z+  Equation  (19)  yields 


C  ^  JS 

5. 


so  that  Equation  .(lU)  for  small  is 

For  leurge  Xz'*',  even  if  z*"  is  small.  Equation  (I8)  becomes 


R(p 


rr  i>  +  c 


Expanding  the  cosines  of  differences  of  angle  occurring  in  Equations  (21)  and 
(22)  yields  the  simultaneous  equation 

Cr*a-0^ 
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(25) 


from  which  D  and  E  are  evaliiated.  Therefore,  Equation  (l8)  is 

Eij») .  Ji?  -V\{  '4^' ;  -  ^  (i?  -  if )  (ii^;  (24) 


As  z+  -►  0  the  product  Ji/3  ^/3  z'*’®/®]  0,  hut  the  product  in¬ 

volving  J  ,  becomes  constant.  Therefore,  the  coefficient  of  J  must  be 

-1/3  -1/3 

zero  if  R  =  0  at  2“*"  =  0.  The  values  of  must  therefore  be  given  by 

%  'ti  •  ♦  ...  (^5^ 

The  equations  for  Rn  aire  therefore 
for  small  r"*"  (center  of  pipe) 

(26) 


for  meditim  r"*" 


R,  (V)  . 


A" 


* 


SAC  aim 


and  for  small  z'*'  =  1  -  r'''  (near  the  wall) 

R.(j*) .  ,.,r 


(27) 


(28) 


Equations  (2lv)  to  (28)  contain  all  the  Infonnation  essential  to  the 
problem  solution.  The  coefficients  in  Equation  (J)  are  found  from  Equation 
(24)  in  accordance  with  Equation  (5)»  Thus  it  is  found  that 


(10  ,  (-0 

y.. 


4**  r(%) 


-n.  .  o,  /, 


and  therefore 


c  ,  (-.r  c’*  *, 


The  derivative  of  R  at  the  wall  (z*»  0)  which  is 


(29) 


(30) 


r(Ej)  3** 


(51^) 


will  be  required  later. 
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Table  I  shows  the  first  ten  eigenvalues  and  the  inportsmt  constants 
for  the  case  of  flow  in  a  round  tube.  Table  II  gives  the  same  data  for  a  flat 
duct  with  opposite  walls  at  the  same  tenperature .  The  development  of  the  flat- 
duct  system  is  similar  to  the  round  duct  and  the  equations  are  given  in  the  ap¬ 
pendix,  numbered  to  correspond  with  the  text. 


TABLE  I  FIRST  TEN  EIGENVALUES  AND  THE  IMPORTANT  CONSTANTS 
FOR  THE  CASE  OF  FLOW  IN  A  ROUND  TUBE 


n 

Cn 

-1/2  C„  R^‘  (1) 

0 

2  2/3 

7.1129 

+1.47989 

0.7303 

1 

6  2/3 

44.489 

-0.80345 

0.53810 

2 

10  2/3 

113.785 

+0.58732 

0.460074 

3 

14  2/3 

215.121 

-0.474993 

0.413743 

4 

18  2/3 

348.457 

+0.404448 

0.381785 

5 

22  2/3 

513.793 

-0.3553^5 

0.357853 

6 

26  2/3 

711.129 

+0.318858 

0.338988 

7 

30  2/3 

940.465 

-0.290488 

0.323555 

8 

34  2/3 

1201.8 

+0.267691 

0.310596 

9 

38  2/3 

1495.1 

-0.248895 

0.29950 

c .  (-,r 

r(\> 

IT 

s  (-0  Ji.f^fCoC  K 

-  A  Vo  _  . 

1  97i  ^ 

a  ’ 

1  V  — 

nr 

Km  • 

»,  J 

1,  ... 

0  m 

?  Cm,  % 

* 

-♦i 
•  T 

z  Sf 

w-t.) 

The  previously  known  eigenvalues  given  by  Jakob  are  shown  in  Table 
III  for  coBopeurison.  Since  the  solution  presented  here  is  valid  for  large 
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m.j' 


TABLE-  II  FIRST  TEN  EIGENVALUES  AND  OHE  IMPORTANT  CONSrANTS  FOR  THE  CASE 

OF  FLOW  IN  A  FLAT  DUCT  WITH  OPPOSITE  WALLS  AT  THE  SAME  TajPERATURE 


n 

Kn 

-Kn  Ifn'  (1) 

0 

1.667 

2.779 

+0.503 

.683 

1 

5.667 

32.11 

-0.121 

.454 

2 

9.667 

93.^5 

+0.0648 

.380 

3 

13.67 

186.9 

-0.0431 

.338 

k 

17.67 

312.2 

+0.0319 

.311 

5 

21.67 

469.6 

-0.0253 

.291 

6 

25.67 

658.9 

+0.0207 

.274 

7 

29.67 

880.3 

-0.0174 

.262 

8 

33.67 

1134 

+0.0150 

.251 

9 

37.67 

i4i9 

-0.0131 

CVl 

• 

K  - 

, 3^ 

•  t-i)  _ 

«(-/;***  at/J 

TT*^ 

yJ(o  =  - 

f .  i"'*  revs') 

0,v» 

TT  3*4 

0,  1,  • 

0  « 

Jlx*)  .  Z  K,  X.'(,) 


TABLE  III  COMPARISON  WITH  PREVIOUSLY  lOICWN  EIGENVALUES 


Results  Obtained 


1  Sellars,  Trihus,  Klein  I 

1  Jakob  { 

1  Analogue 

C(»t>uter 

n 

Cn 

ha. 

Cn 

KmOUIl 

Cn 

0 

2.667 

+1.47989 

0.7303 

2.705 

+1.477 

0.749 

2.71 

1.46 

0.735 

1 

6^667 

-0.80345 

0.5381 

6.66 

-0.810 

0.539 

6.69 

-0.809 

0.533 

2 

10.667 

+0.58732 

0.4601 

10.3 

+0.385 

0.179 

10.62 

+0.592 

0.444 

3 

14.667 

-0.47499 

0.4137 

14.67* 

-0.479* 

14.58 

-0.51 

0.398 

*Attri'buted  to  Lee,  Nelson,  Cherry  and  Boelter 
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and  In  viev  of  the  agreement  even  at  moderate  values  of  Xq,  it  has  been  con¬ 
cluded  that  all  the  eigenvalues  and  functions  are  now  sufficiently  accurately 
known. 

The  heat  flux  at  the  wall  is  computed  from  the  equation 


_V 

( 0  e  ,) 


(52) 


Equation  (32)  is  presented  in  the  above  form  to  bring  it  into  agreement  with 
Jakob . 1 


ARBITRARY  WALL-TEMPERATURE  VARIATIONS 


If  the  wall-tenperature  variation  is  given  by  tw(x),  then,  as  shown 
by  Tribus  and  Klein, 2  the  principle  of  superposition  may  be  applied  and  the 
solution  may  be  written  in  a  Fourier-type  Stieltjes  integral 


t-*. 


(53) 


where  0  is  the  solution  to  Equation  (2)  defined  by  Equation  (3).  The  tenpera- 
ture  of  the  wall  and  fluid  for  x+  <  0  is  to-  The  Stieljes  integral  in  Equation 
(35)  is  evaluated  by  substituting  (dtw/d|)  dS  for  dtw  wherever  tw  is  continuous 
and  substituting  [l  -  9  (x+  -  r+)l[t  (|i+)  -  t  (li")]  as  the  contribution 
of  the  integral  wherever  t^  (x'*’)  has  a  discontinuity.  (See  Tribus  and  KLeinZ 
for  a  more  detailed  discussion.)  The  heat  flux  is  computed  from 


HEAT  FLUX  AT  THE  WALL  GIVEN 


The  inverse  problem;  namely,  "Given  the  heat  flux  at  the  wall,  what 
is  the  temperature?",  may  be  solved  with  the  aid  of  the  Laplace  transform  theory 
Define  the  following  transforms 


(s)  a  -rrs, 0 


(35) 

(36) 


J  frpm  Noni.other.al  Surface.,  by  M.  Tribu.  and 

during  fh-c  ^  A  Sy.po.iu.  held  at  the  Univer.ity  of  Michigan 

,f 
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sx* 

=  5**C  '  -  «•’ 

(37) 

-S** 

•  ”^+  (i^  0  ■  -  0  *■ 

(38) 

„W)  =  5^ 

(39) 

Applying  the  Faltung  theorem  to  Equations  (35)  and  (3^)  yields; 

-r(s.  A*)  *  F(  S,  aO  S  (S^  (*^0) 

(41) 

Q(S)  =  H(53s 

If  the  heat  flux  is  finite,  tw(x+)  will  be  continuous.  Eliminating  tw(8)  from 
the  above  equations, 

^  f(s  aO 

T(S,A*)  ^  ^ 

H(sj  ^ 

and  let  g(x‘'',  r'*')  be  the  inverse  transform  of  G(s,  r"*").  Then,  for  arbitrary 
heat  flux  at  the  wall,  the  tenperature  is  given  by 


Now  define 


(43) 


t- 


(44) 


Thus,  the  problem  is  reduced  to  finding  g(3C^,  r"*"),  which  is  given  by 

.C*  i* 


a(ot*  A*)  *  -L  f  IfL 


M(S) 


AS 


Returning  to  Equations  (37)#  (38)#  and  (3)  it  is  found  that 

C« 


r(s,A*;  *  ±  .  V- 


-  z  ^ 

'»»«•  3+Vj^ 


(45) 


(46) 

(47) 


Because  F  and  H  have  poles  at  s  *  -X%#  the  quotient  F/H  has  no  poles 
except  at  S  s  0  and  the  zeroes  of  H(s)  and  the  zeroes  of  H(s)  must  be  found  nu> 


11 


merically.  Because  H'(s)  is  mouotonic,  it  is  found  that  the  zeroes  of  H(s)  oc¬ 
cur  between  the  Letting  he  the  values  satisfying  H(-7^)  =  0,  from 

the  theory  of  residues 


z 


x:  hv-o 


(*^8) 


Table  IV  gives  the  values  of  7^^^,  H'(-7®jjj)  for  the  first  three  values 
of  m.  The  term  H(0)  has  been  shown  by  others  i  to  be  given  by 


H(o)  s  + 


hence,  the  wall  tenperat\ire  may  be  easily  calculated  with  the  aid  of 


H 


z.  _ 


(*^9) 


(50) 


TABLE  IV  THE  VALUES  OF  7\,  FOR  THE  FIRST  THREE  VALUES  OF  m 

Roots  of  H(8)  =  0,  Values  of  H'(-72jj) 


= 

m 

-1 

7  m 

10  =  -A. 

1 

25.659 

8.854  X  10-3 

4.405 

K  = 

•far.  +  g  0^  2,^  ..  . 

2 

84.624 

2.062  X  10-3 

5.7308 

H'(si 

3 

176.40 

9.435  X  10"* 

6.0084 

* 

/»v  (s-f  Xi  )** 

TABLE  V  VALUES  OF  V2,  H'  (-72)  FOR  THE  FIRST  THREE  VALUES  OF  m  FOR  A  FLAT  DUCT 

Ifl 

Roots  of  H(8)  =  0,  Values  of  H' (-T^jjj) 


H(s)  s  -2  . 

*  Z 


xlt; 

K^XjO) 

7sTvI7 


®  72 

m 


-1 

7®mH’(-72ni) 


1  49. 5^+5  7.1^5  ^  10--* 

2  185. 9**  1.67  X  10-4 

3  1*09.45  6.89  X  10-5 


27.2 

32.1 

35.4 
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A  SAMPLE  CALCULATION  FOR  CONSTANT  WALL  HEAT  FLUX 

By  way  of  illustration  consider  the  computation  of  the  asynptotic  val¬ 
ue  of  the  Nusselt  modulus  for  the  case  of  constant  heat  flux  at  the  wall.  Com¬ 
bining  Equations  (ii-4)  and  (48)  with  q(|)  =  q  =  constant,  the  following  is  ob¬ 
tained. 


t  aV  -t.  -  [f  -z  '  5  -  VJ;  J  (51) 


Letting  Px  »  x"*”,  where  p  =  itk/2WCp,  and  integrating  Equation  (51)  gives 

-t.  .  Ji-  f  -r  -L- 


(52) 


-  Z 


I  -  fc 


-<.\x 


which  may  be  rewritten  as 


=;  j[A.  r  irr  r  —  ‘ 

^  S  vjwv-v) 


i  -  < 


-  XJ  h7-0  "  ? 


r-.] 


(55) 


Equation  (53)  shows  that  far  down  the  pipe  (x"^  -►  «)  the  derivative  of 
t  with  respect  to  x  is  Independent  of  x  or  r'*';  i.e,. 


"  “v 


•f*  p  x'* 


Substituting  this  quantity  into  Equation  (2)  leads  to 

Which  may  be  Integrated  directly  to  give 

t(xlA*)  -  t(X*cJ  ,  ^  ) 

^  f  lU  ^ 

Now  the  mixed  mean  temperature  along  the  pipe  is  given  by 


(54) 


(55) 


(56) 


\*l  i%  '  A  • 


>/  c, 


but  the  mixed  mean  temperature  is  also  defined  by 


(57) 
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=  (  ttf  Cp 


w  ti 


Substituting  Equation  (56)  into  Equation  (58)  and  integrating  results  in 


Combining  Equations  (59)  an*!  (56) 

-t  A*"'  \  .  T  -  3. 

'  I -r  -  -75-  ^  “ 

and  substituting  Equation  (57)  into  Equation  (60) 


at  r+  =  1.  This  expression  reduces  to 


(58) 


(59) 


(60) 


(61) 


(62) 


but  from  Equation  (52),  since  Rn(l)  =  0, 


Hence, 


t,<0 


^  yj  ^  -o.+4-f 


(65) 


(64) 


(Note  that  the  first  three  terms  sum  to  approximately  -0.27<)  Now,  the  Nusselt 
modulus  is  given  by 


Nu.  = 

From  Eqxiation  (60) 

which  when  substituted  into  Equation  (65)  gives 

Hu  » 


(65) 


(66) 


(67) 


Substitution  of  (64),  (58)  and  (55)  into  (65)  yields  the  local  value 
of  the  Nusselt  Modulus  for  the  case  q(x¥)  s  constant. 
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Wu 


(68) 


"  ^  ^  'L  v/h7-^) 


CALCULATION  FOR  LINEARLY  VARYING  WALL  TEMPERATURES 

In  similar  fashion  the  xise  of  the  boimdary  condition  T^ix*)  -  Tq  = 
A  X*  where  A  =  any  constant,  gives: 


7it 


Nl 


a  ■ .  1.  ^ 

(69) 

-t/l  Z 

(70) 

^  ^  K 

e 

(71) 

U  S 


74*  ^  a  -  ^ 


v: 


APPROXIMATIONS  FOR  SMALL  ^ 

Whenever  x"*"  is  small,  a  large  number  of  the  terms  in  the  series,  Equa¬ 
tion  (3)  must  be  taken.  The  Leveque  solution  is  a  good  approximation  for  such 
cases.  As  shown  by  Tribus  and  Klein, ^  the  wall  tenqperature  and  heat  flux  for 
such  a  case  are  related  by 


and 


'■f'o  ' 
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For  flat  ducts 


for  round  ducts 


(ii)  =  ^“7t 


(72) 


(75) 


(74) 


(75) 


Substitution  of  Equation  (75)  into  Equations  (72)  or  (75)  (®uttd  noting 
that  the  mixed  meeui  temperature  of  the  fluid  ds  essentially  equal  to  its  inlet 
value  at  small  values  of  x*)  gives  for  the  three  cases  imder  consideration: 


(76] 


For  constant  wall  temperature; 


V»  rc'M 


♦  -'1 


For  constant  heat  flux 


PK/s) 


X  5  0.00/ 


u  s 


«  /.4373  X 


For  linearly  varying  wall  temperature 


a.oiii  X 


^  ■  S 


(77] 


(78] 


Figure  1  shows  a  graph  of  the  functions  Rq,  R^,  and  Rg  compared  with 
solutions  given  by  Jakob.  Figure  2  shows  the  variations  in  Nusselt  modulus  for 
three  cases 


3 

*Uber  di«  Warmeleitungsfabigkeit  von  Flusnigkeiton, *  by  L.  Graetz. 
Annalen  dor  Phyoik  Cbea.,  vol.  25.  1885,  pp.  337-357. 
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(1)  wall  tesiperature  constant, 

(2)  heat  flvix  constant,  and 

(5)  wall  tenperature  Increasing  lineeurly  along  the  pipe  wall. 

The  Nusselt  modulus  is  defined  hy  the  equation 

N.  •  ^  (79) 

The  mixed -mean  ten?)erature,  twim,  is  determined  by  integrating  the  heat  flux 
from  the  origin  (X+  =  0)  to  the  position  where  q(x)  is  known. 

CONCmSIOWS 

The  methods  used  in  this  paper  have  a  wide  applicability  r  exam¬ 
ple,  the  liquid  metals  systems  analyzed  Poppendiek4  could  be  trt  d  by  the 
methods  used  here.  The  unsymmetrical  booindary  conditions  treated  by  Yih  and 
CermakS  can  also  be  readily  treated  by  these  methods. 

The  authors  are  somewhat  sturprised  at  the  fact  that  whereas  the  as- 
ynptotic  formulae  are  all  supposed  to  be  valid  only  for  very  large  in  ac- 

tmllty  values  of  n  as  small  as  k  seem  to  give  excellent  results .  The  ref  sons 
for  the  good  results  are  not  now  clear. 

APPENDIX  A 


The  equations  for  a  flat  duct  system  with  walls  at  y  =  +  b 


(A-1) 


Defining  Re  =  UUn  pb/^,  x+  =  (xA))(RePr)-S  y*  =  (y/b)  the  equation  to 
be  solved  is 


which  has  a  solution 


1  H 


(A-2) 


mmo  ^ 


(A-3) 


satisfying  ©  «  1  at  x+  «  0,  ©  ♦  0,  r**  -►  0,  if  y(y+)  satisfies 


Part\  • “  Th.r.al  Eatranc.  Region.. 

"*  P®PP*“di«k.  Oak  Ridge  Nationol  Laboratory  Tenn 
ORNL.913.  Phyeic,  eerie.  A.  March,  llsi.  «*»orarory.  xenn.. 

,  r-  ^  Convection  in  Pipe,  ond  Duct..*  by  C.  S.  Yih  ond 

J.  E.  Civil  Engineering  Oepartnent.  Colorado  Agricultural  and 


Y%  (A-**) 

With  Y'(0)  »  Y(l)  =  0,  Y(0)  =  1.  Xq  Is  the  value  of  X  to  permit  Y^^d)  =  0.  The 
coefficients  are  given  by 


x.X 


(A-5) 


?y  the  methods  in  the  text  the  WKB  approximation  is  found  to  be 

"3*" 

for  0  ^  y'*'  <  1. 

Defining  z  =  1  -  y"*",  the  solution  of  A(A-4)  for  z  «  1  is  found  to  be 


(A-14) 


r(}) 


Lie  ^  ~  oi  ;  lox'  2  X  15  iuuna  w  ue 


The  eigenvalues  are 


3‘» '  j'»  H’Q  a 

K.  .  (-0” 


(A-25) 

{A-29) 


(A-30)- 


*  3^  r(%) 

, .  -A  (t.  -t.,  2  e- 

To  obtain  the  fluid  temperature  for  a  given  heat  flux  use 

if 

-- 1 5^’*-s.y^nv*5  (A.4M 


(A-4U) 


The  integrating  kernel,  g,  is  given  by 


X  A 


-<V 


(A-48) 


where  the  -7^  are  the  zeroes  of 

H  (3)  s  -  2;  f  0 


(A-47) 


